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[57] ABSTRACT 

A satellite attitude control system is usable in the ab- 
sence of any inertial yaw attitude reference, such as a 
gyroscope, and in the absence of a pitch bias momen- 
tum. Both the roll-yaw rigid body dynamics and the 
roll-yaw orbit kinematics are modelled. Pitch and roll 
attitude control are conventional. The model receives 
inputs from a roll sensor, and roll and yaw torques from 
reaction wheel monitors. The model produces esti- 
mated yaw which controls the spacecraft yaw attitude. 

8 Claims, 4 Drawing Sheets 
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tion derived solely from the skew gyro is contaminated 
GYROLESS YAW CONTROL SYSTEM FOR A with pitch and roll information. In order to derive yaw 
THREE AXIS STABILIZED, ZERO-MOMENTUM attitude information, from the skew gyro, the roll and 
SPACECRAFT pitch gyros are enabled to provide corrective signals. If, 
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The Government has rights in this invention pursuant 
to the contract NAS5-29600 with NASA. 

BACKGROUND OF THE INVENTION 

This invention relates to attitude control systems 10 
generally, and more particularly, to attitude control 
systems for stabilizing 3-axis zero-momentum satellites 
in the absence of attitude reference information from 
yaw and/or roll inertial references such as gyroscopes. 

Modem satellites are widely used for communica- 15 
tions and Earth sensing applications. All such satellite 
applications require control of the satellite’s orientation 
in space so that sensors and antennas may be pointed in 
appropriate directions. Satellites orbiting the Earth, or 
other heavenly body (hereinafter Earth), do not main- 20 
tain a single face directed toward that body without 
additional control. This control is normally termed 
attitude control. One type of attitude control includes 
the use of one or more momentum wheels, which stabi- 
lize attitude by providing momentum bias or gyroscopic 25 
stiffness. A momentum bias can only indirectly stabilize 
satellite axes lying in a plane orthogonal to the bias axis, 
but direct control is not provided. To satisfy the more 
stringent orientation requirements associated with mod- 
em satellite missions, direct control of the satellite axes 30 
of rotation is required. Three axes are commonly used: 
the yaw axis is oriented toward the Earth, the pitch axis 
is aligned with the satellite’s orbit normal, and the roll 
axis completes the right-handed orthogonal axis set. 
Those skilled in the art know that other non -orthogonal 35 
sets can be used, and that simple transformations relate 
such sets to the orthogonal set. 

The manufacture and launch of satellites is very capi- 
tal intensive. Consequently, to keep the unit cost of 
satellite services low, the satellite must be operated for 40 
a long time. For this reason, satellite reliability is a 
major concern requiring such strong measures as redun- 
dancy, qualification and pre-launch test. 

Typical three axis stabilized satellite attitude control 
is accomplished by directly sensing the three orthogo- 45 
nal attitudes: Yaw(x), Roll(y), and Pitch(z), and com- 
manding corrective control torques through reaction 
wheel actuators, or other torque generators. Such con- 
trol often uses a Earth Sensor Assembly (ESA) to pro- 
vide roll and pitch attitude information. A gyroscope is 50 
used to provide the inertial yaw axis attitude informa- 
tion, because yaw errors (the yaw degree of freedom) 
are not observable with the ESA. Pitch and roll gyros 
are often used in ascent and Earth acquisition phases of 
the satellite’s launch. Once the Earth has been acquired, 55 
however, the roll and pitch attitude information is de- 
rived from the ESA rather than from the gyros. 

The ESA assembly can be fabricated without moving 
parts, and consequently may be very reliable. Even so, 
redundancy assures availability of roll and pitch attitude 60 
information throughout the satellite’s lifespan. How- 
ever the yaw gyro is a mechanical device that is prone 
to failure. Direct redundancy is not normally used be- 
cause of the gyro’s high cost. One additional skew gyro 
provides the only redundancy for the three orthogonal 65 
gyros. 

If the yaw gyro fails, redundant skew gyro could be 
used to provide yaw attitude information, but informa- 


however, one of the roll or pitch gyros (or both) should 
also fail, no useful yaw information is available and the 
satellite attitude undesirably becomes unstable. 

In low cost satellites in which precision pointing is 
not required, the cost of the ESA and gyros may be 
excessive, and other inertial references may not be 
available. An improved attitude control system is de- 
sired. 

SUMMARY OF THE INVENTION 

An attitude control system according to the invention 
controls the orientation of a zero momentum satellite 
relative to a heavenly body. The satellite has gyro- 
scopic freedom in the orbital plane, and derives orienta- 
tion references in the absence of operative gyroscopes. 
First and second reaction wheels have a component of 
momentum along the satellite roll and yaw axes, respec- 
tively. An Earth sensor assembly is mounted on the 
satellite for generating signals representing roll attitude. 
First and second reaction wheel drivers coupled to the 
first and second reaction wheels, respectively, acceler- 
ate the wheels under the control of torque command 
signals representing torque commands about the satel- 
lite yaw and roll axes. First and second tachometers 
coupled to the first and second reaction wheels, respec- 
tively, generate speed signals representing their angular 
velocities. A roll attitude controller is coupled to the 
Earth sensor, and to at least one of the first and second 
reaction wheel drivers for accelerating at least an asso- 
ciated one of the reaction wheels for at least propor- 
tional and rate compensation of the roll attitude. A roll 
reaction wheel torque translator is coupled to at least 
one of the first and second tachometers for differentiat- 
ing the speed signals to generate roll reaction wheel 
torque representative signals. A yaw reaction wheel 
torque translator is coupled to at least one of the first 
and second tachometers for differentiating the speed 
signals to generate yaw reaction wheel torque represen- 
tative signals. A yaw estimator is coupled to the Earth 
sensor, to the roll and yaw reaction wheel torque trans- 
lator and to at least one of the first and second reaction 
wheel drivers, for modeling the linearized roll/yaw 
orbit kinematics of the satellite, and for additionally 
modeling the roll/yaw rigid body dynamics of the satel- 
lite, and for applying signals representing roll attitude, 
and yaw and roll reaction wheel torque representative 
signals to the model, for generating and applying signals 
representing torque commands about the satellite yaw 
axis, to at least one of the first and second reaction 
wheel drivers. 

BRIEF DESCRIPTION OF THE DRAWING 

FIG. 1 is a simplified overall view of an orbiting 
satellite; 

FIG. 2 is a simplified phantom view of a portion of 
the satellite shown in FIG. 1; 

FIG. 3 is a simplified block diagram of a satellite 
attitude control system in accordance with the inven- 
tion; 

FIGS. 4 and 5 are simplified block diagrams of por- 
tions of the system of FIG. 3; 

FIG. 6 is a simplified block diagram of roll and yaw 
rigid body dynamics and orbit kinematics models in 
accordance with the invention; and 
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FIG. 7 is a plot of the yaw performance of a satellite 
using a gyroless yaw control system in accordance with 
FIG. 7. 

DESCRIPTION OF THE INVENTION 5 

In FIG. 1, a satellite 8 includes a body 10. Satellite 8 
orbits about the Earth or another heavenly body 12 
Body 10 supports first and second solar panels 14 a and 
146. Body 10 supports an Earth sensor 16 which is 
capable of sensing the horizon, thereby providing pitch 10 
and roll attitude signals. Earth sensor 16 is oriented 
along an axis 18 which is parallel to body yaw axis 20. 
Body roll and pitch axes 22 and 24, respectively, are 
mutually orthogonal and orthogonal to yaw axis 20. 
The orbital direction of motion of satellite 8 is along 15 
body roll axis 22, and pitch axis 24 is orthogonal to the 
orbital plane. 

In FIG. 2 a pitch reaction wheel 218 has its momen- 
tum axis parallel to (in this case coincident with) pitch 
axis 24. A pitch wheel driving assembly illustrated as a 20 
block 220 a controls wheel 218 in a conventional man- 
ner. 

Also in FIG. 2, a yaw reaction wheel 214 with its 
momentum axis parallel to (in this case coincident with) 
yaw axis 20, and is driven by a driving assembly or 25 
motor 216c. The speed of the yaw wheel 214, is moni- 
tored by a tachometer 2166. A roll reaction wheel 210 
with its momentum axis parallel to (in this case coinci- 
dent with) roll axis 22, is driven by a driving assembly 
or motor 212c. The speed of the yaw wheel 214, is 30 
monitored by a tachometer 2126. 

In general, the gyroless yaw control is provided by a 
model of the satellite roll-yaw rigid body dynamics, as 
in the prior art, together with a further model of the 
small angle roll-yaw orbit kinematics. The body dynam- 35 
ics are the rotational accelerations of the satellite body 
about its own center of mass in response to externally 
applied torques, which torques may be generated by an 
internal apparatus such as a reaction wheel. Small angle 
orbit kinematics describe the satellite’s orbital motion. 40 
While it would be desirable to use general orbit kine- 
matics in the model, the complexity of the non-linear 
calculations is not justified in view of the relatively 
small accuracy increment. 

FIG. 3 is a simplified block diagram of a zero- 45 
momentum 3 axis stabilized satellite attitude control 
system using rigid body dynamic and orbit kinematic 
models in accordance with the invention. Elements of 
FIG. 3 corresponding to those of FIGS. 1 and 2 are 
designated by the same reference numerals. In FIG. 3 50 
Earth sensor 16, produces measured or actual pitch and 
roll signals on data paths 350 and 352 respectively for 
application to pitch and roll control compensators 330 
and 332 respectively. Details of pitch control compen- 
sator 330 are described in conjunction with FIG. 4 and 55 
those of roll control compensator 332 are described in 
conjunction with FIG. 5. Pitch and roll compensators 
330 and 332 produce pitch and roll torque commands 
respectively, which are applied to pitch and roll wheel 
drivers 220a and 212a, respectively. Pitch and roll 60 
wheel drivers 220a and 212a drive pitch and roll reac- 
tion wheels 218 and 210, respectively. The angular 
speeds of pitch and roll reaction wheels 218 and 210 are 
measured by pitch and roll tachometers 2206 and 2126, 
respectively. 65 

A system model 310 in FIG. 3 includes a model, 
illustrated as a block 312, of the roll-yaw rigid body 
dynamics of the particular satellite. Such a model in- 
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eludes information relating to at least the inertia about 
the roll and yaw body axes. System model 310 also 
includes a model of the small angle roll-yaw kinematics, 
illustrated as a block 314. This information is not spe- 
cific to any characteristic, such as mass or inertia, of the 
satellite itself. Rather, only the satellite orbit rate oq is 
significant. Inputs to system model 310 include the ac- 
tual roll derived from Earth sensor 16, and applied over 
data path 352, and roll and yaw torques, which are 
derived from roll tachometer 2126, and from an yaw 
tachometer 2166, by way of roll and yaw translators 
illustrated as blocks 316 and 318 respectively. This mini- 
mal information allows the yaw to be estimated. The 
estimated yaw signal produced by system model 310 is 
applied to a yaw control compensator 334, similar to 
pitch and roll compensators 330 and 332, respectively, 
to produce yaw torque commands which are applied to 
a yaw reaction wheel 214, by way of a yaw wheel 
driver, 216a. Thus, the yaw attitude may be corrected 
without any direct measurement thereof. 

Prior art arrangements also include yaw estimators, 
as described for example in U.S. Pat. No. 4,521,855, 
issued Jun. 4, 1985 in the name of Lehner, et. al. While 
such systems estimate yaw, their models couple the roll 
and yaw dynamics by the pitch bias momentum. Thus, 
they can only be used for satellites including a pitch 
momentum wheel providing momentum bias. Under 
many conditions, the pitch bias momentum results in 
coupling between the roll and yaw attitude motion. 
Such coupled motions complicate the control, and may 
limit achievable yaw and roll attitudes. For example, to 
change the yaw orientation in a pitch bias momentum 
satellite, application of a yaw torque results in the need 
to apply a compensating roll torque due to the gyro- 
scopic coupling. Ordinarily, this undesirably results in 
expenditure of additional fuel. 

In FIG. 4, the measured or actual pitch (0), is applied 
from Earth sensor 1 16 (FIG. 3) over data path 350, in 
parallel, to a scalar multiplier 410 for multiplication by 
a factor Kpz, to an integrator 414, for integration and 
multiplication by a factor Kiz, and to a differentiator, 
416, for differentiation and multiplication by a factor 
Kdz. The resulting proportional, integral, derivative 
(PID) signals are summed in a summing circuit 412. 
While PID control is preferred, PD control has been 
found to be satisfactory for many control applications in 
which steady state sensor errors are not important. The 
pitch PID signal at the output of summer 412 is ampli- 
fied, and if needed, analog-to digital-converted, in an 
amplifier, 220a. The amplified signal is applied to pitch 
wheel 218. In a particular embodiment of the invention, 
further described below, the pitch wheel dynamics are 
described in the Laplace domain by the transfer func- 
tion: 75 s/(600 s+ 1), meaning that the amplified signal 
voltage results in a proportional pitch body torque for 
amplified signal frequencies greater than 1/600 sec. 

In FIG. 5, the measured or actual roll (<f>), is applied 
from Earth sensor 16 (FIG. 3) over data path 352, in 
parallel, to a scalar multiplier 510 for multiplication by 
a factor Kpy, to an integrator 514, for integration and 
multiplication by a factor Kiy, and to a differentiator, 
516, for differentiation and multiplication by a factor 
Kdy. The resulting proportional, integral, derivative 
(PID) signals are summed in a summing circuit 512. The 
roll PID signal at the output of summer 412 is amplified, 
in an amplifier 212a. The amplified signal is applied to 
roll wheel 210. In the abovementioned embodiment of 
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the invention, the roll and pitch reaction wheel dynam- 
ics are identical. 

In the more detailed block diagram of FIG. 6, the 
actual roll signal (<J>) from Earth sensor 16 of FIG. 3 is 
applied over data path 352, to a non-inverting (+) input 
port of a summer, 610. Summer 610 also receives an 
estimated roll signal, derived as described below, at its 
inverting (— ) input port. The output of summer 610 is 
the system error, termed “estimator” error. The estima- 
tor error is applied over a data path 612 to scalar multi- 
pliers 614, 616, 644, and 646, for multiplication by se- 
lected constants. The estimator error scaled by multi- 
plier 614 is applied to a 4- input port of a summer 618 
together with actual yaw torques as translated from the 
yaw reaction wheel of FIG. 3 by yaw translator 318, 
and the roll gyroscopic torque produced by the scalar 
multiplier 648. Similarly, the estimator error scaled by 
multiplier 616 is applied to a -f input port of a summer 
620 together with actual roll torques as translated from 
the roll reaction wheel of FIG. 3 by roll translator 316. 
The yaw gyroscopic torque produced by the scalar 
multiplier 650 is applied to a (— ) input port of 620. The 
output signals from summers 618 and 620 are the total 
yaw and total roll body torques, respectively. 

The rigid body dynamics of the satellite, designated 
generally as 312 in FIG. 6, and more specifically desig- 
nated 622, are represented by blocks 624, 626, 628, 630, 
648 and 650. The total roll body torque, produced at the 
output of summer 620 is applied to a scalar multiplier 
626 for dividing the total roll body torque by the known 
satellite roll axis rotary inertia, (Iy) to produce roll 
acceleration, a) y . Similarly, the total yaw body torque 
produced at the output of summer 618 is applied to a 
scalar multiplier 624, for dividing the total yaw body 
torque by the known satellite yaw axis rotary inertia, 
(lx) to produce yaw acceleration, ai*. The yaw and roll 
body accelerations produced at the outputs of multipli- 
ers 624 and 626, respectively, are applied to time inte- 
grators 628 and 630, respectively, to produce yaw and 
roll body angular velocities co*, and & y respectively. 
The output of the roll-yaw rigid body dynamics 312 are 
the estimated roll and yaw body angular velocities a>* 
and a)y respectively. The yaw body angular rate pro- 
duced at the output of integrator 628 is applied to a 
scalar multiplier 650 for multiplying the yaw body an- 
gular rate by the product of the satellite orbital rate (107 
o) and the difference between the known satellite pitch 
rotary inertia (I 2 ) and roll inertia (!>,) to produce gyro- 
scopic roll torque. The roll body angular rate produced 
at the output of integrator 630 is applied to a scalar 
multiplier 648 for multiplying the total roll body angu- 
lar rate by the product of the satellite orbital rate (a>o) 
and the difference between the known satellite pitch 
rotary inertia (I 2 ) and yaw inertia (I*) to produce gyro- 
scopic yaw torque. Block 648 represents the coupling of 
rotation about the roll axis into yaw torque attributable 
to the orbital rotation of the satellite about the Earth. 
Block 650 similarly represents the coupling of rotation 
about the yaw axis into roll torque attributable to the 
orbital rotation of the satellite about the Earth. 

The roll-yaw orbit kinematics estimator 314 in FIG. 6 
includes a yaw summing circuit 632, which receives the 
estimated yaw body angular velocity o) x , at a + input 
port, together with estimator error scaled by multiplier 
644. Similarly, estimator 314 of FIG. 6 includes a roll 
summing circuit 634 which receives the estimated roll 
body angular velocity a>y, at a 4- input port, together 
with estimator error scaled by multiplier 646. Summer 
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632 produces the estimated yaw attitude rate i|f, which 
differs from the body yaw rate, o) x in that the yaw atti- 
tude rate includes the effect of yaw-roll orbital cou- 
pling. Similarly, summer 634 produces the estimated 
roll attitude rate <f>, which differs from the body roll 
rate, a^in that the roll attitude rate includes the effect of 
yaw-roll orbital coupling. 

The outputs of summers 632 and 634 are applied to 
time integrators 638 and 636, respectively, to to inte- 
grate the estimated yaw and roll attitude rates, respec- 
tively, to produce estimated yaw and roll attitudes, iff 
and respectively. 

The estimated roll at the output of time integrator 636 
is fed back to a -f input of summer 632 by way of a 
scalar multiplier 640 which multiplies by a factor o>o, 
representative of the satellite’s orbital angular velocity. 
The estimated roll produced at the output of integrator 
636 is applied to an inverting input of summer 610. 
Similarly, the estimated yaw at the output of time inte- 
grator 638 is multiplied by too in a scalar multiplier 642 
and applied to an inverting (— ) input of summer 634. 
The estimated yaw produced at the output of integrator 
638 is the desired output of system model 310. 

The estimated yaw output from system model 310 is 
applied to yaw control compensator 334 of FIG. 3, 
which is identical to the pitch or roll compensators 330 
or 332, respectively. 

In order to understand the concepts underlying the 
invention, one may consider that the estimator is used to 
control the yaw attitude. The essence of the estimator 
can be described by a simplified, four degree of freedom 
mathematical model of the roll/yaw orbit kinematics, 
and the roll/yaw rigid body dynamics of an orbiting 
vehicle. The degrees of freedom can be written as a 
state vector, X, as follows: 

JT=[Yaw (Iff), Roll (<J>), Yaw Rate (a> x ) Roll Rate 
(<•>>)] 

The differential equations for the four states are: 

ip^<Ojc+oiod 

<f>=o) y —o)o iff 

U x =h~ l [T x +(I z -Iy)<0OMy] 

U>y=Iy- l [T X ~(Ig-I x )o> 0O) X ] (1) 

Where I**, lyy, and I^are the vehicle yaw, roll and pitch 
rotary inertias about its center of mass, T*and T^are the 
total yaw and roll torques, and o>o is the orbital angular 
velocity. 

Since the state equations are cross-coupled, it is math- 
ematically possible to reconstruct the estimated yaw 
attitude, and body rate, cu*, by knowing only the 
actual satellite roll attitude, <f>. This permits using the 
four degree of freedom model as an estimator of the 
actual satellite yaw, roll, yaw rate and roll rate states, or 
degrees of freedom. The estimator forms the difference 
between the actual satellite roll (sensed by its Earth 
sensor or by other means) and the estimated roll to 
produce the estimator error. Then it drives the model 
with the estimator error to converge its state toward the 
actual satellite state. The convergence is accomplished 
by adding the estimator error, scaled by estimator gains 
(Ke(l) through Ke(4)) to each of the four differential 
equations of motion for the four states. 
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An actual embodiment of the invention was imple- 
mented on an operating satellite with a 100 minute orbi- 
tal period, and having identical yaw, roll and pitch 
reaction wheels, the dynamics of which are defined 
above, and also having identically tuned yaw, roll and 5 
pitch attitude control compensators (330, 332, and 334 
of FIG. 3) adjusted to 0.1 rad/sec bandwidth, and in 
which the the estimator gains were —960.0, 0.609, 
-0.384, and -0.899 for Ke(l), Ke(2), Ke(3), and Ke(4) 
respectively. When the gyroscopically actuated yaw 10 
control system was disabled, and the above described 
estimator was given control of the satellite yaw attitude, 
the yaw attitude error as a function of time was plotted. 
FIG. 7 is a copy of the telemetered data. In FIG. 7, the 
plot designated generally as 710 represents degrees of 15 
yaw attitude error. In FIG. 7 control was transferred to 
the inventive yaw estimator at about 72 min along the 
time scale. Starting from zero error, the system con- 
verges and stabilizes the yaw attitude at about 80 min- 
utes (portion 712 of plot 710). Region 714 represents the 20 
estimator’s control in response to unmodelled environ- 
mental disturbance torques. Interruption 716 in plot 710 
represents a data drop-out. At about 98 minutes, the 
estimator was re-started with zero initial conditions as a 
test. Again, for reasons that are not fully understood, 25 
the estimator experiences a brief transient in the oppo- 
site direction. Following 1 10 minutes, the system again 
stabilizes, but was disturbed in the region 720, by the 
satellite’s automatic magnetic momentum unloading 
system. Such unloading takes place periodically during 30 
normal satellite operation because of the accumulation 
of disturbance momentum. In the interval between 120 
and 150 minutes, the unperturbed estimator yaw control 
can be seen as a region 722 of plot 710. Following 150 
minutes, another automatic magnetic momentum un- 35 
loading event occurs which obscures the estimator re- 
sponse. It appears from the unperturbed portions of the 
plot that the yaw error does not exceed about three 
degrees when operating using estimated yaw control. 
The constant negative offset of about 1.5 degrees is 40 
attributed to an unmodelled constant aerodynamic body 
fixed yaw torque. 

Other embodiments of the invention will be apparent 
to those skilled in the art. For example, analog and 
digitally operated controls may be freely interchanged 45 
as conditions require. Also, the estimated yaw rate used 
in the yaw control compensator 334 of FIG. 6 can be 
derived by differentiating the estimated yaw as de- 
scribed, or it may instead be taken directly from the 
output of integrator 628. The estimator gains may be 50 
calculated by using a numerical pole-placement algo- 
rithm, by using a Linear Quadratic Regulator algo- 
rithm, or by any other method. 

What is claimed is: 

1. An attitude control system for controlling the ori- 55 
entation of a zero momentum satellite relative to a heav- 
enly body, said satellite having gyroscopic freedom in 
the orbital plane, and deriving orientation references in 
the absence of operative gyroscopes, comprising: 

a satellite body defining roll and yaw axes, said yaw 60 
axis being Earth pointing, said roll axis being paral- 
lel to the satellite motion, said satellite body having 
yaw and roll body inertias about said yaw and roll 
axes respectively; 

first and second reaction wheels mounted on said 65 
body, each of said first and second reaction wheels 
having a component of momentum along the satel- 
lite roll and yaw axes, respectively; 
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Earth sensing means mounted on said body for gener- 
ating signals representing roll attitude of said satel- 
lite relative to said heavenly body; 

first and second reaction wheel driving means cou- 
pled to said first and second reaction wheels, re- 
spectively, for accelerating said wheels under the 
control of torque command signals representing 
torque commands about said satellite yaw and roll 
axes; 

first and second tachometer means coupled to said 
first and second reaction wheels, respectively, for 
generating speed signals representing the angular 
velocities of said reaction wheels; 

roll attitude control means coupled to said Earth 
sensing means, and to at least one of said first and 
second reaction wheel driving means for process- 
ing said roll attitude error signals with at least pro- 
portional and rate compensation of said roll atti- 
tude, for generating said torque command signals 
representing said torque commands about said roll 
axis, for accelerating at least an associated one of 
said reaction wheels; 

yaw attitude control means coupled to at least one of 
said first and second reaction wheel driving means 
for processing estimated yaw signals with at least 
proportional and rate compensation of said yaw 
attitude, for generating said torque command sig- 
nals representing said torque commands about said 
yaw axis, for accelerating at least an associated one 
of said reaction wheels; 

roll reaction wheel torque translation means coupled 
to at least one of said first and second tachometer 
means for differentiating said speed signals to gen- 
erate roll reaction wheel torque representative 
signals; 

yaw reaction wheel torque translation means coupled 
to at least one of said first and second tachometer 
means for differentiating said speed signals to gen- 
erate yaw reaction wheel torque representative 
signals; 

yaw estimating means coupled to said Earth sensing 
means, to said roll and yaw reaction wheel torque 
translation means, and to at least one of said first 
and second reaction wheel driving means, for mod- 
elling the linearized roll/yaw orbit kinematics of 
said satellite, and for additionally modelling the 
roll/yaw rigid body dynamics of said satellite, and 
for applying said signals representing roll attitude, 
and said yaw and roll reaction wheel torque repre- 
sentative signals to said model, for generating and 
applying signals representing torque commands 
about said satellite yaw axis, to said at least one of 
said first and second reaction wheel driving means. 

2. A system according to claim 1 wherein said reac- 
tion wheel axes are non-orthogonal, and further com- 
prising: 

a momentum-axis transformation means coupled to 
said tachometer means and to said wheel drive 
means for transforming coordinates between said 
non-orthogonal axes to orthogonal axes. 

3. A system according to claim 1 wherein said first 
reaction wheel is a yaw reaction wheel, and said one of 
said first and second reaction wheels to which said yaw 
estimating means is coupled is said yaw reaction wheel. 

4. A system according to claim 1 wherein said satel- 
lite further defines a pitch inertia about a pitch axis, and 
said yaw estimating means includes roll/yaw rigid body 
dynamic modelling means comprising: 
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yaw and roll estimating means for generating an esti- 
mated yaw and roll signals in response to estimated 
roll and yaw body rate signals and scaled estimated 
roll error signals; 

first differencing means coupled to said Earth sensing 5 
means, and to said yaw and roll estimating means 
for generating said estimated roll error signals, 
representative of the difference between said roll 
representative signal and estimated roll signals; 

first and second scaling means coupled to said first 1 ° 
differencing means for multiplying said estimated 
roll error signals by first and second scalar quanti- 
ties, respectively, to produce first and second 
scaled estimated roll error signals. 

first summing means coupled to said first scaling 15 
means and to said yaw reaction wheel torque trans- 
lation means for summing said yaw reaction wheel 
torque representative signals, said scaled estimated 
roll error signals, and an estimated yaw gyroscopic 
torque signal, for generating estimated yaw torque 20 
signals; 

second summing means coupled to said second scal- 
ing means, and to said roll reaction wheel torque 
translation means, for summing said roll reaction ^ 
wheel torque representative signals, said scaled 2 
estimated roll error signals, and an estimated roll 
gyroscopic torque signal, for generating estimated 
roll torque signals; 

yaw body rate estimating means coupled to said first 
summing means, for generating estimated yaw 
body rate by integrating the quotient of said esti- 
mated yaw torque signals divided by said satellite 
yaw inertia; 

roll body rate estimating means coupled to said sec- 35 
ond summing means, for generating estimated roll 
body rate by integrating the quotient of said esti- 
mated roll torque signals divided by said satellite 
roll inertia; 

yaw-roll cross-coupling torque estimating means cou- 4 ^ 
pled to said roll body rate estimating means and to 
said first summing means, for multiplying, by the 
orbital rate, the difference between said pitch and 
roll body inertias, to produce estimated yaw gyro- 
scopic torque signals, and for coupling said esti- 45 
mated yaw gyroscopic torque signals to said first 
summing means; 

roll-yaw cross-coupling torque estimating means cou- 
pled to said yaw body rate estimating means and to 
said second summing means, for multiplying, by 50 
the orbital rate, the difference between said pitch 
and yaw body inertias, to produce estimated roll 
gyroscopic torque signals, and for coupling said 
estimated roll gyroscopic torque signals to said 
second summing means. 55 

5. A system according to claim 4 wherein said yaw 
estimating means includes roll/yaw orbit kinematics 
modelling means comprising: 

third and fourth scaling means coupled to first differ- 
encing means for producing third and fourth scaled 60 
estimated roll error signals; 

third summing means coupled to said roll/yaw rigid 
body dynamic modelling means and to said third 
scaling means for summing together said estimated 
yaw signals, said third scaled estimated roll error 65 
signals, and a yaw attitude kinematic correction 
signal to produce estimated yaw attitude rate sig- 
nals; 
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fourth summing means coupled to said roll/yaw rigid 
body dynamic modelling means and to said fourth 
scaling means for summing together said estimated 
roll signals, said fourth scaled estimated roll error 
signals, and a roll attitude kinematic correction 
signal to produce estimated roll attitude rate sig- 
nals; 

third integration means coupled to said fourth sum- 
ming means, and to said first differencing means, 
for integrating said estimated roll attitude rate sig- 
nals, to produce said estimated roll attitude signals; 

fourth integration means coupled to said third sum- 
ming means, and to said first differencing means, 
and to said yaw attitude control means for integrat- 
ing said estimated yaw attitude rate signals, to pro- 
duce said estimated yaw attitude signals; 

roll and yaw orbital rate scaling means coupled to 
said third and fourth integrating means, respec- 
tively, and to said third and fourth summing means, 
respectively, for multiplying said estimated roll and 
yaw respectively by signals representative of said 
orbit rate for producing said yaw and roll attitude 
kinematic correction signals, respectively. 

6 . A method for controlling the orientation of a zero 
momentum satellite relative to a heavenly body, said 
satellite having gyroscopic freedom in the orbital plane, 
and deriving orientation references in the absence of 
operative gyroscopes, the satellite body defining roll 
and yaw axes, said yaw axis being Earth pointing, said 
roll axis being parallel to the satellite motion, said satel- 
lite body having yaw and roll body inertias about said 
yaw and roll axes respectively, said method comprising 
the steps of: 

spinning first and second reaction wheels mounted on 
said body, to produce components of momentum 
along said satellite roll and yaw axes, respectively; 

generating signals representing roll attitude of said 
satellite relative to said heavenly body; 

generating torque command signals representing said 
roll attitude about said roll axis; 

applying said torque command signals to said at least 
one of said first and second reaction wheels; 

generating roll reaction wheel speed signals repre- 
senting the effective angular velocities of said reac- 
tion wheel; 

generating yaw reaction speed signals representing 
the effective angular velocities of said reaction 
wheel; 

differentiating said roll reaction wheel speed signals 
to generate roll reaction wheel torque representa- 
tive signals; 

differentiating said yaw reaction wheel speed signals 
to generate yaw reaction wheel torque representa- 
tive signals; 

modelling the linearized roll/yaw orbit kinematics of 
said satellite for producing a roll/yaw orbit kine- 
matics model; 

modelling the roll/yaw rigid body dynamics of said 
satellite for producing a roll/yaw rigid body dy- 
namics model; 

applying said signals representing roll attitude error, 
and said yaw and roll reaction wheel torque repre- 
sentative signals to said roll/yaw rigid body dy- 
namics model for producing estimated roll and yaw 
body rate signals; 

applying said roll and yaw body rate signals, and said 
roll error signals to said to roll/yaw orbit kinemat- 
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ics model for generating estimated yaw and roll 
attitude error signals; 

generating yaw torque command signals representing 
said estimated yaw error about said yaw axis; 
applying said yaw torque command signals to at least 
an associated one of said reaction wheels, for accel- 
erating said wheel, for generating yaw torque. 

7. A system according to claim 1 wherein said satel- 
lite further defines a pitch inertia about a pitch axis, and 
wherein said step of modelling said roll/yaw rigid body 
dynamics of said satellite comprises the further steps of: 
taking the difference between said roll representative 
signal and estimated roll signals to produce signals 
representing estimated roll error; 
multiplying said estimated roll error signals by first 
and second scalar quantities, respectively, to pro- 
duce first and second scaled estimated roll error 
signals. 

summing said yaw reaction wheel torque representa- 
tive signals, said scaled estimated roll error signals, 
and an estimated yaw gyroscopic torque signal, for 
generating estimated yaw torque signals; 
summing said roll reaction wheel torque representa- 
tive signals, said scaled estimated roll error signals, 
and an estimated roll gyroscopic torque signal, for 
generating estimated roll torque signals; 
generating estimated yaw body rate by integrating 
the quotient of said estimated yaw torque signals 
divided by said satellite yaw inertia; 
generating estimated roll body rate by integrating the 
quotient of said estimated roll torque signals di- 
vided by said satellite roll inertia; 
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multiplying said estimated roll body rate, by the 
product of orbit rate and the difference between 
said pitch and roll body inertias, to produce esti- 
mated yaw gyroscopic torque signals; 

5 coupling said estimated yaw gyroscopic torque sig- 
nals to said first summing means; 
multiplying said estimated yaw body rate, by the 
product of orbit rate and the difference between 
said pitch and yaw body inertias, to produce esti- 
10 mated roll gyroscopic torque signals; 

coupling said estimated roll gyroscopic torque signals 
to said second summing means. 

8. A method according to claim 7 wherein said step of 
modelling said roll/yaw rigid body dynamics of said 
15 satellite comprises the further steps of 

scaling said estimated roll error signals to produce 
third and fourth scaled estimated roll error signals; 
summing together said estimated yaw signals, said 
third scaled estimated roll error signals, and a yaw 
20 attitude kinematic correction signal to produce 
estimated yaw attitude rate signals; 
summing together said estimated roll signals, said 
fourth scaled estimated roll error signals, and a roll 
attitude kinematic correction signal to produce 
25 estimated roll attitude rate signals; 

integrating said estimated roll attitude rate signals, to 
produce said estimated roll attitude signals; 
integrating said estimated yaw attitude rate signals, to 
produce said estimated yaw attitude signals; 

30 multiplying said estimated roll and yaw signals, re- 
spectively, by signals representative of said orbit 
rate for producing said yaw and roll attitude kine- 
matic correction signals, respectively. 
***** 
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